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TITLE 

TYPE 

NUMBER OF W0RD6 

TEMPORARY STORAGE 

PRESET PARAMETERS 



ACCURACY 
DURATION 



DESCRIPTION 



ENTRY 



Linear Ifatrix Equation Solver and General Matrix Inversion 
Closed routine with two auxiliaries ^ D01 or S* 1 * 01 ) 

115 

0, 1, 2, 3, b t 5, 6, 7 (also see below) 

S3, S4. During input of routine the following parameters 

are required: 

Memory positions 3 = 00 F 00 q^^ F 

k = 00 F 00 qg F 

where q. and q^ are the first locations of auxiliary 
routines I and II respectively. 

A function of the order and conditioning of the matrices, 
n. = order of matrix A, m = number of columns of matrix B 

- 11 + kn - 3n 2 + Oln 3 (m = n) n 3 /500 (m = 1) 

Upon appropriate entry, this routine solves the linear 
matrix equation AX = B or inverts the matrix A, where 
the matrices A and B must satisfy the conditions: 

a) A is non-singular and of size (n x n) 

b) B is of size (n x m) 

The magnitudes of n and m are governed by the available 
storage space. In addition to the routine itself, 

[ n j n " t " f + n + m + nm] memory positions are required 

directly after the routine. Also, consideration must 

be given for storing the original matrices A and B and 

finally the solution X. (See note on Auxiliaries). 

A distinction is to be made between the case when B = I 

(i.e. when ve wish to invert A) and the case when B f I, 

in which we are in essence solving m sets of n simultaneous 

equations in n unknowns . 
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1) B =.I: 



Entry; 



p 


JO n 
L5 P 


p + 1 


26 x 
00 n 



2) B £ I: Entry; 



p 


50 n 
L5 P 


p + 1 


26 x 

00 m 



AUXILIARIES 



where s n = order of A 

m = columns of B 

x = first location of this routine. 

The auxiliaries provide a flexible means for handling 



the matrices A and B and solution X. By this method one can: 

1) successively generate A and B row "by row, so that 
the whole matrix need not he in the machine at one 
time. 

2 ) have the matrices located on the drum, Williams 
memory, or on tape. 

One can, in the same way, obtain successive columns 
of X and store them on the drum or punch them on tape so that the complete 
matrix X need not necessarily be stored in its entirety at any one time. 

Auxiliary I: Exit to this subroutine is made from 
word 22 of the main routine. This auxiliary must provide successive rows of 
A and B as follows: 

a) B = I. The identity matrix is automatically 
generated in the main routine so that only 
rows of A must be entered, starting at memory 
position (x + 115) to (x + llU + n). 

b) B f I. Successive rows of the augmented matrix 
[A,Bj must be provided and read into memory 
positions starting at (x f 115) to (x + 11^ + n + m) 
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Counting of number of rows (n) is not necetsaary in the 
auxiliary nor must provision be made for a link since direct transfer is to 
be im&e at the end of the auxiliary to the right hand side of memory position 

(x + 21). 

The matrices must be scaled so that no element is > l/2. 

Auxiliary II. Since each column of X is computed 
independently, this auxiliary provides opportunity for disposing of that column 
before the next one is obtained. Exit to this subroutine is made from word 
(x + 105) j the elements of the respective columns and scaling factor being found 
in memory position (x + 115 ) to (x +115 + n). No counting of number of 
columns (m) is necessary in the auxiliary, nor must provision toe made for a link. 
Direct transfer is made at the end of the auxiliary to the left hand side of 
memory position (x ■♦ ■ 106). 

NOTE Following the n elements of each column «?? * is a scaling 

factor which may be different for each column* It locates the position of the 
decimal point as lying after the position occupied by the number J . For example, 

the sequence 

+ 2000 - 05C0 4- 2500 4- 1000 

- ^000 + OUy) - 2kOQ + 0100 

■+■ 0500 + 0100 + 0050 -* iOuO 

represents the third order xr-atrl* 

2 -30 -3 

- .5 .? .1 

2,5 ^ ,05 

-.1 
In case we are inverting A the solut'lo;-. a • must 

n+1 ' 

be multiplied by 10 , if the original matrix A was scaled by 10 . 

.MATHEMATICAL METHOD See Routine M 15 fcr mathematical method. The ' precautions 
noted there for detecting sinularities of A are incorporated in this routine. 
Instead of printing out F it will cause a machine hangup franc memory position 
(x 4 1.1.1) where order ON 3276 M 3277 is .located. 



DAT S 5/31/55 s KT: 5/29/56 
PROGRAMMED BY W. L. Frank 
APPROVED BY 




Igr 



LOCATIOI 


ORDER 


10X03 PiOE 1 





00 K(ML*) 

42 15L 




Set p 




40 F 




Start 


1 


41 3F 








41 4F 




•* Clear counters 


2 


41 7F 
F5 F 






3 


42 10L 




Set p + 1 




L4 87L 




Right address p + 2 


4 


46 j6h 
46 83L 






5 


42 108L 




Link 




10 20F 




Set address of y = x + 114 + a 


6 


42 75L 
1*2 91L 






7 


42 19L 
LO 43L 






8 


42 3F 
L5 45L 




Store n 


9 


42 22L 








L5 91L 




Set y 


10 


42 16L 








L5 (p+l)F 


3 


. 


11 


42 4F 




Set m 




L4. 91L 






12 


42 23L 
42 45L 






13 


42 112L 

42 114L 




- Set address y + m - t 


11* 


00 20F 








46 113L 






15 


41 6F 




Clear counter 




L5 (p)F 





Test for layers ion or solving Az = B 


16 


36 21L 








41 (y)F 


10,17* 





M Ifc 




M l -k 



| LOCATION 


ORDER 


NOTES PAGE 4 


51 


L5 XF 


^ 


•■ 


below the diagonal 




40 (t)F 


26,52 r 


c f A 


52 


F5 51L 




1 


■ 


te 5X1* 






53 


42 ^5L 






■ 


L5 **8L 




! 


54 


L4 109L 




* T 




40 48L 




1 

B t 


55 


k6 50L 




1 ^ 




42 44l 


■ 




56 


LO 115L 




! 


57 


32 kkL 
13 5F 




If ze?d ctonH descale 


"■' 


32 6% 




' 


58 


IX 5F 




>i/n 




32 6.U. 






59 


L5 66l 








k6 6% 




' 


60 


& 6ki 

r?6% 




J If 'N(5) > 1/2 multiply row by l/?. 
T (scale 4pwa) 


61 


42 65L 




t 
I 




1*5 5F 




! 


• 6a, 


00 IF 

1 




*" Determine If possible to scale up 




ko 5? 




<l/2 N(5),<l/4 


63 


LL 5F 








32 60L 






64 


50 7P 
L5 (x)F 


68,28' 






65 


10 (1)F 


34*,59 l 






■• , ■ 


00 (1)F 


35,61 




* , „ 


66 


50 2F 




Waste 






40 (*)F 


V* £ jjife^r 




Rescale rows 


67 


F5 66lr~" 
42 66L 


t 




- 

" 


68 


te 64L 








; ' ' ■ 


LO U4L J 









*t Ik 



LOCATIOI 



ORDER 



69 
TO 
71 
72 
73 



75 

76 

17 
78 
79 

ao. 

81 
32 
83 

85 



36 64l 
F5 6f 
4o 6f 
L5 7F 

LO 6F 
36 2^L 

F5 7F 

40 7F 
LO 5F 
32 7^L 
22 8L 

41 5F 
L5 11JJ* 
^0 (y)F 

L3 (y)F 
36 111L 

41 6? 
kl 7F 
L5 ^5L 
k2 88L 
LO 4F 
IA 5F 

42 831, 
L5 88L 
LO kF 
k2 88L 
L5 83L 
46 88L 
50 (y)F 

71 ( )? 
40 F 
L5 83L 
FO 7F 
LO kf 



SOTES 
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M \h 



4' 
80,86 



J 



Determine If row I must have further 
eliminations (done i times) 



- Count nuaaher of rows 

Repeat for next row 
Set counter 

Set scaling factor 



Terminate calculation if scaling < 2 



Clear counters 



o-39 



i 1 
j-4 



- Set addresses 



Loounor 

86 

87 
88 

89 

90 
91 
92- 
93 
9* 
95 
96 

97 

98 

99 
100 
101 

10B 



ORDER 

k2 85L 
22 92L 
85 llfcL 
50 IF 
50 (y)F 
7* ( )F 
Lfc F 
1*0 F 
LL F 
32 92L 
50 111L 

7J (y)F 

22 75L 
L5 8qL 
LO 109L 
kO 88L 
te 98L 
te 101L 
h6 103L 
F5 6F 
kO 6> 
L5 7F 
LO 6F 
36 87L 
kl 6> 
L3 (a lt )F 
36 111L 

L6F 

36 91L 

26 101L 

I$F 

66 (a tl )F 

22 102L 

81 F 



soras 



BUS 6 



Calculate E 
J-i+1 



Vj 



M i4 



Also acta as constant 



82^93 



6' 



S* 



9* 1 



> 1/2? 



Rescale and start again 



Sat addresses 



Count (n»l+l)tiines for row i 



Reset counter 

Snd If zero on diagonal 



Test If division Is proper 
Waste 



Waste 



IOCAKON 
103 

ICft 

105 
106 

107 
108 

109 
110 
111 

113 

Ilk 



ORDER 

kO (y.-.l)7 95 

F5 7F 

^0 TF 

LO 3F 

36 S4 

22 80L 

F5 5F 

40 5? 

LO k¥ 

32 108L 

26 75L 

26 (pf2)F 

00 IP 

00 IP 

7L ^095F 

LL ^095F 

00 F 

00 1000 0000 0000 J 

S6 21L 



NOTES 
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M .\4 



kl (t)P 
NO (t)F 
L3 F 
JO 2F 
40 (t)F 



13 
1* 



13 1 



. Count n rms 

Exit to Auxiliary i: 
Repeat 

■ Count m cdltnazus 

Repeat 

Link 

1-2 7J 
1/10 
.Sod constants 



H 



—j 



